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ABSTRACT 

Emission from X-ray binaries (XRBs) is a major component of the total X-ray luminosity of normal galaxies, 
so X-ray studies of high-redshift galaxies allow us to probe the formation and evolution of XRBs on very long 
timescales (~10 Gyr). In this paper, we present results from large-scale population synthesis models of binary 
populations in galaxies from z = 0 to ~20. We use as input into our modeling the Millennium II Cosmological 
Simulation and the updated semi-analytic galaxy catalog by Guo et al. to self-consistently account for the star 
formation history (SFH) and metallicity evolution of each galaxy. We run a grid of 192 models, varying all the 
parameters known from previous studies to affect the evolution of XRBs. We use our models and observationally 
derived prescriptions for hot gas emission to create theoretical galaxy X-ray luminosity functions (XLFs) for several 
redshift bins. Models with low common envelope efficiencies, a 50% twins mass ratio distribution, a steeper initial 
mass function exponent, and high stellar wind mass-loss rates best match observational results from Tzanavaris & 
Georgantopoulos, though they significantly underproduce bright early-type and very bright (L x > 10 41 ) late-type 
galaxies. These discrepancies are likely caused by uncertainties in hot gas emission and SFHs, active galactic 
nucleus contamination, and a lack of dynamically formed low-mass XRBs. In our highest likelihood models, we 
find that hot gas emission dominates the emission for most bright galaxies. We also find that the evolution of the 
normal galaxy X-ray luminosity density out to z = 4 is driven largely by XRBs in galaxies with X-ray luminosities 
between 10 4() and 10 41 erg s _1 . 

Key words: binaries: close - galaxies: stellar content - stars: evolution - X-rays: binaries - X-rays: diffuse 
background - X-rays: galaxies 
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1. INTRODUCTION 

X-ray binaries (XRBs) are believed to be major contributors 

to the overall X-ray luminosity of normal galaxies (those not 

dominated by the emission of a nuclear supermassive black hole; 
BH; Fabbiano 1989, 2006; Kim & Fabbiano 2003). Normal 

early-type galaxies have older stellar populations and their 
X-ray emission is dominated by low-mass XRBs (LMXBs) and 
hot interstellar medium (ISM). On the other hand, the X-ray 

emission of normal late-type galaxies, which are still actively 
forming stars, have significant contributions from both LMXBs 

and high-mass XRBs (HMXBs). 

X-ray and multiwavelength studies of galaxies using Chandra 

and XMM-Newton have yielded a great deal of information 
about the X-ray luminosities of galaxies, including many X-ray 
correlations that have been established to hold out to at least 
z — 1 (e.g., Lehmer et al. 2007, 2008; Symeonidis et al. 2011; 
Vattakunnel et al. 2012; Basu-Zych et al. 2013; Cowie et al. 
2012). These relations include a strong correlation between 
X-ray emission from HMXBs and the star formation rate (SFR) 
of galaxies (e.g., Ranalli et al. 2003; Gilfanov et al. 2004; 
Lehmer et al. 2010; Mineo et al. 2012) as well as a scaling 
relation between the emission from LMXBs and the stellar 
mass of a galaxy (Gilfanov et al. 2004; Gilfanov 2004; Lehmer 


et al. 2010; Boroson et al. 2011; Zhang et al. 2012). Recent 
ultradeep Chandra and multiwavelength surveys (e.g., Brandt 
& Hasinger 2005) have allowed for robust tests of these relations 
in very distant galaxies. Lor example, Basu-Zych et al. (2013) 
use a 4 Ms exposure of CDF South (Xue et al. 2011) and X-ray 
stacking to study faint X-ray sources out to z ~ 8, finding that 
the relation between X-ray production and SFR undergoes a 
small amount of evolution out to z ~ 4 that is likely driven by 
the metallicity evolution of HMXBs. 

Galaxy X-ray luminosity functions (XLFs) derived from 
recent observations show significant evolution with red- 
shift (Norman et al. 2004; Ranalli et al. 2005; Ptak et al. 
2007; Tzanavaris & Georgantopoulos 2008). Tzanavaris & 
Georgantopoulos (2008, hereafter T&G08) use data from three 
Chandra deep fields (CDF-S, E-CDF-S, and CDL-N) and the 
wide-area survey XBootes to compile observations of 207 X-ray 
luminous normal galaxies (101 early-type and 106 late-type) 
out to z ~ 1.4. They find a clear evolution of the galaxy XLL 
normalization with redshift that is driven almost exclusively 
by late-type galaxies. More specifically, this evolution is pro- 
portional to (1 + z) k with k — 2.2 ± 0.3 for the total popula- 
tion, k = 2.4 4l 2Q for late-type galaxies, and, for the early-type 
population, k — 0.7*\g (consistent with zero). Because XRBs 
are major contributors to the total X-ray emission of normal 
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galaxies, observationally derived XLFs can put constraints on 
theoretical models of XRB formation and evolution. 

At present, there has been little theoretical work done on 
the evolution of XRB populations over cosmological timescales 
(White & Ghosh 1998; Ghosh & White 2001; Zuo & Li 201 1). 
It is thought that XRBs could play a major role in the evolution 
of these XLFs (White & Ghosh 1998; Ghosh & White 2001; 
Norman et al. 2004). Ghosh & White (2001), using a semi- 
empirical, semi-analytical approach, linked XRB lifetimes with 
SFRs, showing that SFRs that are evolving on cosmological 
timescales significantly affect the XRB populations and, there- 
fore, the integrated galactic X-ray emission. This predicted evo- 
lution should be evident even at lower redshifts (z < 1 ; White 
& Ghosh 1998). 

Recently, the advances in available multiwavelength obser- 
vations of distant galaxies, as well as our understanding of bi- 
nary stellar evolution and galaxy formation and evolution, have 
reached a level of maturity that allows us to conduct an in-depth 
study of the XRB populations of distant galaxies. In this paper, 
we study the evolution of XRBs on cosmologically significant 
timescales, using data from detailed, large-scale simulations. 
We use data from a catalog created by Guo et al. (2011) using 
semi-analytical galaxy evolution models applied to the recent 
Millennium Cosmological Simulation. These data are used in 
tandem with the binary population synthesis (PS) code, Star- 
Track, to simulate the XRB populations of individual galaxies 
from z = 0 to ~20, taking into account the full star formation 
and merger histories of each galaxy. From these models, we de- 
rive the integrated X-ray emission of each galaxy and compare 
the resulting galaxy XLFs and their evolution to observed galaxy 
XLFs. Our goal is to obtain better constraints on the parameter 
space for models of XRB formation and evolution, and to better 
understand the nature of the X-ray emission of galaxies at high 
redshifts. 

Recently, Fragos et al. (2013, hereafter F13) used similar 
techniques as those described in this paper to study the evolution 
of the global XRB population with redshift. They model how the 
total universal specific X-ray luminosities (Lx per unit stellar 
mass, SFR, and volume) of LMXBs and HMXBs evolve over 
cosmic time out to z ~ 20. Their models were constrained by 
observed luminosities of HMXB and LMXB populations in the 
local universe. They found that the LMXB population dominates 
the total population at low redshifts, with HMXB contributions 
becoming dominant for redshifts higher than z ~ 3.1. 

The outline of the paper is as follows. Section 2 describes our 
simulation tools, StcirTrack and The Millennium Simulation II, 
and the methodology we follow in developing our models of 
XRB populations in galaxies. Section 3 describes how we 
compare our models to observational results, namely those of 
T&G08, and the statistical analysis we use to determine our best 
models. In Section 4, we describe and discuss our results, and 
we conclude with a summary in Section 5. 

2. SIMULATING X-RAY LUMINOSITIES OF GALAXIES 
2.1. The Millennium Cosmological Simulation 

The Millennium Cosmological Simulation is an unprece- 

dented computational effort to simulate the dark matter dis- 

tribution in the universe (see Springel et al. 2005 for details). 

In this study, we use the data from the most recent Millen- 
nium Run II (hereafter MRII; Boylan-Kolchin et al. 2009). This 
is an /V-body simulation that follows the evolution of 2160 2 3 * * 
particles, each of mass 6.9 x 10 6 ft" 1 M 0 within a comoving 


box with sides each of size 100 ft" 1 Mpc. The cosmological 
model used in the simulation is a A cold dark matter model 
with Q m = 0.25, Qa = 0.75, Cl/, = 0.045, and h = 
Z/o/100 km s" 1 Mpc" 1 = 0.73. 

The MRII has 60 snapshots in time that were saved and 
analysis was done to identify substructures within the dark mat- 
ter distributions, including dark matter halos. Guo et al. (2011, 
hereafter G1 1) use a semi-analytic procedure to track the evolu- 
tion of the galaxies that exist within these halos. Once subhalos 
are identified, their merger trees are derived. The evolution of 
these subhalos provides the base for the galaxy formation model. 
The models used by G1 1 build upon the work of De Lucia & 
Blaizot (2007), making improvements to the treatment of su- 
pernova (SN) feedback, reincorporation of ejected gas, galaxy 
sizes, the distinction between satellite and central galaxies, and 
the effect of the environment on galaxies. While semi-analytical 
models do not supply accurate details about individual galax- 
ies, they are very useful for understanding general character- 
istics of large populations of galaxies. These semi-analytical 
models, when applied to the MRII simulation, are able to accu- 
rately reproduce observed characteristics of galaxy populations, 
e.g., the abundance and large-scale clustering of low-z galaxies, 
the Tully-Fisher relation, stellar mass and luminosity functions 
of low-z galaxies, the halo-galaxy-mass relationship, and the 
evolution of the cosmic star formation density. However, these 
models overproduce passive low-mass galaxies and fail to repro- 
duce the observed abundances, clustering, and mass functions 
of high-redshift (z > 1.0) galaxies. In this paper, we will be 
comparing with X-ray observations of galaxies out to redshift 
1.4, which is in a regime where the Gil model is still fairly 
accurate. 

The result of Gil’s semi-analytic model is a catalog of the 
galaxy population at 60 different times between z ~ 20 (about 
13.4 billion years ago) and the present day. These catalogs 
include properties such as metallicity, stellar mass, bulge mass, 
the mass of hot and cold gas, rest-frame luminosity magnitudes, 
etc., for each of the galaxies in the simulation box, as a function 
of time. 

2.2. StarTrack 

To simulate the XRB populations of the galaxies from MRII, 
we use StarTrack, a current binary PS code that has been 
tested and calibrated using detailed binary star calculations and 
incorporates all the most important physical processes of binary 
evolution (Belczynski et al. 2002, 2008). 

1. The evolution of single stars and non-interacting binary 
components, from their birth, taken as the time of their 
initial emergence onto the main sequence, to compact 
remnant formation using evolutionary formulae of Hurley 
et al. (2000) modified as described in Belczynski et al. 
(2008). Various wind mass-loss rates and their effect on 
stellar evolution are also incorporated into the code and 
have been recently updated (Belczynski et al. 2010). 

2. The time evolution of orbital properties. StarTrack numeri- 
cally integrates a set of four differential equations describ- 
ing the evolution of orbital separation, eccentricity, and 
component spins, taking into account tidal interactions, 
magnetic braking, gravitational radiation, and stellar wind 
mass losses. 

3. All types of mass-transfer phases. This includes both stable 
and unstable mass-transfer processes, which are driven 
by either nuclear evolution or angular momentum loss. 
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Parameter 

Notation 

Value 

Reference 

Initial orbital period distribution 3 

F(P) 

Flat in logP 

Abt (1983) 

Initial eccentricity distribution 3 

F(e ) 

Thermal F(e) ~ e 

Heggie (1975) 

Binary fraction 3 
Magnetic braking 3 

fb in 

50% 

Ivanova & Taam (2003) 

Metallicity 3 

z 

0.0001, 0.0002, 0.005, 0.001, 
0.002,0.005,0.01,0.02, 0.03 


IMF (slope) 3 


-2.35 or -2.7 

Kroupa (2001) and Kroupa & Weidner (2003) 

Initial mass ratio distribution 3 

F(q) 

Flat, twin, or 50% flat plus 50% twin 

Kobulnicky & Fryer (2007) and Pinsonneault & Stanek (2006) 

CE efficiency b 

«CE 

0.1, 0.2, 0.3, or 0.5 

Podsiadlowski et al. (2003) 

Stellar wind strength 3 

*7wind 

0.25, 1.0, or 2.0 

Belczynski et al. (2010) 

CE during HG b 


Yes or No 

Belczynski et al. (2007) 

SN kick for ECS/AIC C NS 3 


20% of normal NS kicks 

Linden et al. (2009) 

SN kick for direct collapse BH b 


Yes or No 

Fragos et al. (2010) 


Notes. 

a Observationally constrained parameters. 
b “Free” parameters. 

c Electron capture supernova/accretion-induced collapse. 


Unstable mass transfer is encountered most often as a 
direct consequence of rapid stellar expansion during nuclear 
evolution, but angular momentum loss may also lead to 
instability. 

4. SN explosions, which are treated by taking into account 
mass/angular momentum losses as well as SN asymmetries 
(through natal kicks to neutron stars (NS) and BHs). 

5. X-ray emission, which is tracked for accreting binaries 
with compact object primaries (both for wind-fed and 
Roche-lobe overflowing systems). The resulting X-ray 
luminosities are calculated from the secular averaged mass 
accretion rate, but are not calculated for unstable accretion 
phases because the timescales are very short. 

The models in this paper include a recent revision of StarTrack 
that incorporates updated stellar winds and their recalibrated 
dependence on metallicity (Belczynski et al. 2010). However, 
two more recent upgrades have not been incorporated into these 
results, as the simulations were run long before the changes 
were made. The first update includes a revised NS and BH 
mass distribution based on fully consistent SN simulations 
(Belczynski et al. 2012; Fryer et al. 2012). The second, most 
recent upgrade improves upon the treatment of donor stars in 
common envelope (CE) events via usage of the actual value of 
the X parameter, the measure of the central concentration of the 
donor and envelope binding energy, for which usually a constant 
value is assumed (Dominik et al. 2012). 

Table 1 lists the input parameters of our PS models, which 
can be put into two categories. In the first category there are the 
parameters that correspond to the initial properties of the binary 
population. These values are relatively well constrained by the 
most recent observational surveys. Also in the group are stellar 
wind prescriptions and natal kick distributions, which can also 
be constrained by observations. In the second group there are 
the truly “free” parameters that correspond to poorly understood 
physical processes, which we are not able to model in detail. One 
of these truly “free” parameters is the CE efficiency (o!ceX which 
measures how efficiently orbital energy loss is transformed into 
thermal energy that will expel the donor’s envelope during the 
CE phase. We note that in our calculations we combine oice 
and X, the binding energy parameter described above, into one 
CE parameter. Whenever we mention the CE efficiency a Cf ., 


we refer in practice to the product ckce x X, effectively treating 
o'er, x a as a free parameter (see Belczynski et al. 2008 for 
details). 

We create a grid of 192 PS models, a subset of those used 
in FI 3, run for 9 different metallicities and each simulating 
5.12 x 10 6 stars for 14 Gyr. In this grid, we varied all the 
parameters known from previous studies to affect the evolution 
of XRBs and the formation of compact objects in general 
(Belczynski et al. 2007, 2010; Fragos et al. 2008, 2010; Linden 
et al. 2009). Specifically, we vary the CE efficiency, initial binary 
mass ratio distribution, initial mass function (IMF), SN kicks for 
direct collapse (DC) BHs, and stellar wind strength. We also take 
into account the possibility of CE inspirals with Hertzsprung gap 
(HG) donors that could terminate binary evolution barring the 
subsequent XRB formation (Belczynski et al. 2007). 

For all models, we assume a Maxwellian distribution of SN 
kicks given by Hobbs et al. (2005), with a = 265 km s -1 . For 
compact objects formed with partial mass fallback, the natal 
kicks given by the Hobbs et al. (2005) distribution are decreased 
by a factor of (1 — //*), where fp, is the fraction of the stellar 
envelope that falls back after the SN explosion. In our standard 
prescription, DC BHs, BHs formed with ffi, = 1, are given 
no natal kick. However, due to recent theoretical evidence that 
even the most massive stellar BHs have probably received small 
asymmetric kicks (Linden et al. 2010; Valsecchi et al. 2010), in 
some models used in this work we set a lower limit (0.1) on the 
amount by which the natal kicks may be decreased due to mass 
fallback, allowing for small natal kicks to be given to DC BHs. 

The mass of the primary star in each binary is determined by 
the adopted IMF. It is important to note here that, because we 
sample the IMF with only the primary star, we are only sampling 
the high-mass end of the IMF because the primary stars that 
form XRBs must be massive enough to form a BH or NS. The 
mass of the secondary star is calculated using a distribution 
function for the binary mass ratio, q — M secon dary/(H pr j mar y. We 
vary the distribution of q between a flat distribution in the range 
q — 0-1 and a distribution that has 50% of the binaries follow 
a distribution with q — 0-1 and the other half follow a “twins” 
distribution, with q — 0.9-1. 

The models and their numbers used in this work are the 
same as those used in F13, except here we exclude the models 
97-192, which have all systems following the pure “twins” 
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^-distribution. FI 3 show that models with the pure “twins” 
distribution are unphysical, as they fail to reproduce local 
populations of XRBs. Thus, we exclude them here. 

We note that our PS code calculates the bolometric luminosity 
of each XRB based on the rate of mass transfer. In order to 
compare our model results with observed data sets we need to 
estimate the X-ray luminosity in a specific energy band, which 
in this study is the soft X-ray band of Chandra (0.5-2.0keV). 
In order to calculate the bolometric correction, we used two 
sets of published X-ray spectra from Galactic NS and BH 
XRBs at different spectral states (Remillard & McClintock 
2006; Wu et al. 2010). Following the same procedure outlined in 
FI 3, we derive to the bolometric correction factors for different 
types of XRBs and use these results to estimate the 0. 5-2.0 keV 
X-ray luminosity of our modeled XRB population. 

We also note that the PS models used here take into account 
only binary systems formed in the field and not those formed 
via dynamical interactions in dense stellar systems. Dynamical 
formation in globular clusters (GCs), for example, is a signif- 
icant formation pathway for LMXBs in old, GC-rich elliptical 
galaxies (e.g., Humphrey & Buote 2008; Zhang et al. 2012). 
Further, the LMXB populations formed in GCs can be as much 
as 2-3 times more luminous than the field population in bright 
elliptical galaxies (Irwin 2005). 

For each model and each metallicity value, we calculate the 
X-ray luminosity as a function of age for single bursts of star 
formation, where we also take into account the effect of transient 
XRBs (see Fragos et al. 2008, 2009). Taking into account the 
assumed IMF and initial mass ratio distribution, we normalize 
the total X-ray luminosity to a nominal population of 10 10 M 0 . 
This quantity, Lx, S pec(0( er g s — 1 / 10 10 M 0 ), is the specific X-ray 
luminosity of a single-age stellar population as a function of its 
age. The specific X-ray luminosity coming from our PS models 
can be convolved with the star formation history (SFH) and 
metallicity evolution of a galaxy to calculate the total X-ray 
luminosity of its complex stellar population. 

2.3. Convolving StarTrack with the Gil Catalog 

The MRII catalog created by G1 1 corresponds to 60 snapshots 
that span a redshift range from z = 0 to ~20. For each galaxy, 
we can derive its complete progenitor tree. Each progenitor 
galaxy has a unique SFR and metallicity, so for every stellar 
population in a target galaxy we know during what time frame 
and at what metallicity that population was created. We then 
convolve the SFHs with Lx, S pec(0 derived from the PS models 
for the appropriate metallicity values. 

The SFRs given for each galaxy in the Gil catalog are 
averaged over the entire time step. At, between subsequent 
snapshots so that the total new stellar mass created in a given 
progenitor galaxy is M new = SFR prog Af. In order to account 
for the possibility of starbursts, we assume that all new stellar 
mass forms in a 20 Myr burst occurring at a random time 
between subsequent snapshots, t, < r burst < f i+1 , where t, is the 
timestamp associated with the snapshot of a given progenitor 
galaxy. This effect is important only for the HMXBs of young 
populations, and the effect on LMXBs is minimal since their 
evolution occurs on timescales much longer than the time steps 
between snapshots. A 20 Myr duration is reasonable, given that 
the most cited values for starburst durations are around 3-10 Myr 
(e.g., Thornley et al. 2000; Harris et al. 2004), though there is 
evidence for longer burst durations on the order of a few 10 8 yr 
(e.g., McQuinn et al. 2010). 


By summing the soft-band X-ray luminosities of all the 
stellar populations in a given galaxy, we derive the integrated 
X-ray luminosity from XRBs. The end result is a catalog 
of integrated XRB luminosities of galaxies within the MRII 
comoving volume from z = 0 to z ~ 20. 

3. COMPARING WITH OBSERVATIONS 
3.1. Galaxy Classification 

When comparing our results to the observations of T&G08, 
we want to distinguish between early- and late-type galaxies. 
For their sample, T&G08 cross-correlate with other surveys 
to obtain optical counterparts for their X-ray-selected galaxies, 
which they used for classification. 

The classification of a galaxy as early or late type can be 
based either on its morphology or its spectroscopic properties. 
In many observational surveys of distant galaxies (e.g., GOODS; 
Giavalisco et al. 2004), the morphologies of most galaxies that 
are observed cannot be determined due to inadequate spatial 
resolution, and therefore colors are used instead. 

For color classification of our model galaxies as early or late 
type, we adopt the method developed by Bell et al. (2004). 
They showed that it is possible to define the population of 
early-type galaxies empirically by using the bimodality of the 
color distribution, which they studied out to z ~ 1 . The MRII 
database includes absolute rest-frame magnitudes in the Sloan 
Digital Sky Survey ugr filters, which can easily be transformed 
to the UBV filters used in Bell et al. (2004). The magnitudes 
include the effects of dust extinction. Following the Bell et al. 
(2004) prescription, we define early-type galaxies to be galaxies 
where (U — V) ^ 1.15 — 0.31z — 0.08 * ( M v — 5 log 10 h + 20). 
Figure 1 shows plots of ( U — V) versus My — 5 log 10 h for 
the MRII galaxies at various redshifts with the cutoff function 
overlaid on top and the bimodality is clearly present. 

If we instead define galaxy types based on morphology with 
late-type galaxies having (M\, u \ ge / M tota i) < 0.7 and early-type 
galaxies (M\, u \ ge /M to t a i) > 0.7, we find that there is approxi- 
mately a 1% and 0.5% contamination among the color defined 
late and early types, respectively. Thus, these methods give 
nearly identical results, but using colors to define morphol- 
ogy allows us to better simulate observations. The morphology 
method of classification, since it is independent of color, pro- 
vides a check on our color classification. The fact that they 
both yield similar results is encouraging and indicates that the 
colors provided by the Gil catalog are able to yield accurate 
morphology classifications. 

3.2. Creating Model X-Ray Luminosity Functions 

For our analysis, we only select galaxies with stellar masses 
greater than 10 s solar masses, as galaxies with mass less 
than this are very unlikely to have X-ray luminosities that are 
observationally relevant. For instance, the dwarf galaxies in the 
SINGS sample, with masses ~10 7 M Q , generally have X-ray 
luminosities below 10 37 erg s _1 (0.5-8 keV) with many having 
no binaries detected at all above 5 x 10 36 erg s -1 (0.5-8 keV; 
L. Jenkins et al. 2013, in preparation). The G1 1 catalog contains 
several hundreds of these more massive galaxies at very high 
redshift (z. = 19.9) and over 2 million at z = 0. It should also be 
noted that we assume that all of the galaxies in our sample are 
normal galaxies. This is a valid assumption because galaxies 
with bright active galactic nucleus (AGN) only constitute 
~2%-4% of all galaxies and, therefore, any selection effects on 
our data would be minimal (Xue et al. 20 1 0; Haggard et al. 20 1 0; 
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-22 -20 -18 -16 -22 -20 -18 -16 -22 -20 -18 -16 
M v - 5 log, 0 h M v - 5 log 10 h M u - 5 log 10 h 

Figure 1. Rest-frame U — V color of simulated galaxies against the absolute magnitude in V-band, My — 5 log 10 h. It is clear that the G1 1 galaxy catalog exhibits the 
same bimodality as observed galaxies. The red line corresponds to 1.15 — 0.3 lz — 0.08 * (My — 5 log 10 h + 20). This may be compared to the red sequence fitting of 
U — V colors by Bell et al. (2004, see their Figure 1). Galaxies that lie above the red line are considered to be early-type galaxies and those that lie below are late-type 
galaxies. 

(A color version of this figure is available in the online journal.) 


Silverman et al. 2009). Lower luminosity AGN have been found 
in a much higher percentage (~30% — 40%) of LINER galaxies 
(Ho et al. 1997). However, since LINER galaxies themselves 
make up only one-fifth to one-third of all galaxies (Ho et al. 
1997), this effect is also rather minimal. Additionally, since this 
is for lower luminosity AGN, it is likely that our more luminous 
galaxies would still be dominated by hot gas and XRB emission 
and would be classified as a normal galaxy. For example, Flohic 
et al. (2006) find that AGN in their sample of LINER galaxies 
contribute only 60% of the 0.5-10 keV luminosity when only 
considering the central regions of the galaxies. 

To compare our results with observations, we derive the XLF 
for our galaxies by calculating the number density of galaxies 
versus their integrated X-ray luminosity. The simulation data 
are all within a single comoving volume of constant size. For 
the time slice represented in each snapshot, <p(L ) is defined as 


0(L) = 


V (L m in, L max ) 


Vmrii<5z, 


( 1 ) 


Here, L mm , L max are the bin limits, Vmrii is the volume of the 
MRII simulation, which is (100 Mpc /i -1 ) 3 , and b/ is the size of 
the luminosity bin in log space, i.e., S L = log 10 (L max /L m i n ). 
Herein lies a major difference between the theoretical and 
observational luminosity functions. Observational surveys study 
a range of redshifts within a light cone. An entire volume of 


space cannot be observed at a constant redshift so a range of 
redshifts is explored. Thus, when calculating <p(L), observers 
such as T&G08 use methods like the one found in Page & 
Carrera (2000) which uses the following definition: 


<P(L, z) 


N 

r, Lm “ r^^dzdL 

" L.min ^ Znan(i-,) dz 


(2) 


Here, dV /dz represents the rate of change of the survey volume 
with respect to redshift and z max (L), z mm (L) are the redshift 
ranges for a source as a function of luminosity such that it 
stays within the flux limits of the survey and within the redshift 
interval. Our simulated galaxies, on the other hand, exist within 
a comoving volume and our snapshots capture all galaxies 
that exist at a given redshift. In order to compensate for this 
difference, we adopt similar redshift intervals as used by T&G08 
and calculate the XLF for each redshift individually using 
Equation (1). This gives us (j){L. z). Then, for each luminosity 
bin centered at L, we take the average value of <p(L, z) over all 
z in the interval, giving us an estimate for cp(L ) for that bin in 
that redshift interval. 


3.3. X-Ray Luminosity from Hot Gas 

In addition to XRBs, the hot ISM in a galaxy can have a 
significant contribution to its overall X-ray luminosity. T&G08 
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Table 2 

Parameters and Likelihood Values for Models Referred to in This Paper 


Model 

«CE a 

IMF Exponent 

^wind* 3 

CE-HG C 

^-distribution* 1 

DC BH kick e 

Rank f 

log(L(0|M)/L ref )S 

205 

0.1 

2.7 

2.0 

No 

50-50 

0.0 

1 

0.0000000 

229 

0.1 

2.7 

2.0 

Yes 

50-50 

0.0 

2 

-0.057250977 

277 

0.1 

2.7 

2.0 

Yes 

50-50 

0.1 

3 

-0.53945923 

245 

0.1 

2.7 

1.0 

No 

50-50 

0.1 

4 

-1.6570358 

253 

0.1 

2.7 

2.0 

No 

50-50 

0.1 

5 

-1.7356873 

273 

0.1 

2.35 

2.0 

Yes 

50-50 

0.1 

6 

-2.3401947 

269 

0.1 

2.7 

1.0 

Yes 

50-50 

0.1 

12 

-5.6114807 

249 

0.1 

2.35 

2.0 

No 

50-50 

0.1 

10 

-4.1068573 

248 

0.5 

2.7 

1.0 

No 

50-50 

0.1 

55 

-47.292496 

197 

0.1 

2.7 

1.0 

No 

50-50 

0.0 

50 

-42.623398 

241 

0.1 

2.35 

1.0 

No 

50-50 

0.1 

59 

-51.771774 

261 

0.1 

2.7 

0.25 

No 

50-50 

0.1 

81 

-92.800415 

53 

0.1 

2.7 

1.0 

No 

Flat 

0.1 

22 

-15.239967 


Notes. 

a CE efficiency parameter. 
b Stellar wind strength parameter. 

c 0: CE from Hertzsprung gap donor allowed and 1 : not allowed. 
d Binary mass ratio distribution. 

e SN kicks given to direct collapse black holes. 0.0 = no SN kick given and 0.1 = small SN kick given. 
f The rank of the model based on the likelihood value. 

g The log of the ratio of the likelihood of the given model to that of the highest likelihood model. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding 
its form and content.) 


do not distinguish between emission from XRBs and hot gas in 
their analysis of X-ray bright galaxies. Their analysis is done 
in the soft X-ray band, so emission from the hot ISM becomes 
important and needs to be taken into account when calculating 
the total integrated luminosities of our galaxies. 

We use observationally derived relations to estimate the 
X-ray luminosity of hot gas in early-type galaxies from their 
X-band luminosity. It has been shown that there is a power- 
law relationship between the X-ray luminosity of the hot 
ISM in early-type galaxies and both the X-band luminosity 
(Lx oc L l k ' m ) of the galaxy and the temperature of its hot 
gas (T oc Lx 214 ; Boroson et al. 2011). The X-band luminosities 
of MRII early-type galaxies are easily calculated from mass 
and age using synthetic stellar population models (Bertelli et al. 
2008). With these relations, we estimate both the full-band X-ray 
luminosity and the temperature of the hot gas in each early-type 
galaxy. The spectrum of hot diffuse gas is assumed to be that 
of a collisionally ionized diffuse gas as calculated by the APEC 
XSPEC model and the ATOMDB code (Foster et al. 2012). The 
gas temperature estimate from the empirical relations is used as 
input to the APEC model in order to calculate the luminosity of 
the hot gas in the soft X-ray band. 

For late-type galaxies, we estimate the hot gas X-ray luminos- 
ity based on the SFRs given in the G1 1 catalog and the power-law 
relationship catalog between the soft-band (0.5-2 keV) X-ray 
luminosity of the hot ISM and the SFR for late-type galaxies 
(L x a SFR 107 ; Strickland et al. 2004a, 2004b). 

In summary, we estimate the total X-ray emission from 
hot gas in all of the galaxies in the Gil catalog. We add 
those values to each galaxy’s integrated X-ray emission from 
XRBs, calculated using StarTrack, to obtain the total integrated 
X-ray luminosity of each galaxy. We find that on average XRB 
emission contributes to ~50%-60% of the 0.5-2 keV emission 
from bright (L x > 10 38 ) late-type galaxies and ~40% of the 
0.5-2 keV emission from bright early-type galaxies for our best- 
fitting model (205). Hence, we find that hot gas emission has an 


appreciable effect on the galaxy XLFs. See Section 4 for more 
details on our results. 


3.4. Statistical Analysis 

From deep Chandra survey observations, T&G08 present 
early- and late-type galaxy counts in several luminosity bins 
and two redshift intervals, 0 ^ z < 0.4 and 0.4 ^ z < 1.4. 
They also provide total galaxy counts split into three redshift 
intervals, but for our analysis we will focus on the early- and 
late-type counts. Associated with each count is a survey volume 
that depends on both luminosity and redshift. Let the set of 
counts be d = {d,j\i = 1, . . . N, j = 1, . . . M } and associated 
volumes be Vjj, where i ranges over the N — 5 luminosity bins 
and j ranges over the M = 2 redshift bins. T&G08 assume that 
the di j are subject to Poisson statistical errors. 

Similarly, for a particular choice of parameters, 9, (for 
example, see Table 2) our model produces a set of counts n = 
{n, j\i = I .... /V , / = 1, . . . M } for galaxies in each luminosity 
and redshift bin in our (lOOMpc /? -1 ) 3 model volume. We 
assume that the counts n are drawn from a Poisson distribution 
with (unknown) means X = {X{j\i = 1, . . . N, j = 1, . . . M } 
(note that this X is separate from the one used before to describe 
the CE efficiency parameter). Because we only observe one 
particular set of counts, n, we do not measure the rates X implied 
by our model directly, but instead must treat X as a nuisance 
parameter whose distribution under the observed n must be 
integrated over. 

Bayes’ rule relates the posterior probability of model 
parameters 9, p(9\d), to the likelihood of the data under the 
model, p(d\9), the prior probability of the model parameters 
before the data have been observed, p(9), and a normalizing 
constant, p(d), called the evidence, that is independent of 9 via 


P(9\d) = 


P(d\9)p(9) 

P(d) 


( 3 ) 
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Models 1-96 Models 193-288 

Figure 2. Likelihood values for all of the models used in this work. L re f is the highest likelihood value among our 192 models. The model numbers correspond to the 
same models used in F13, though here we exclude models 97-192, as F13 found them to be inconsistent with observations. 


Writing the likelihood in terms of the (unknown) true mean X We choose a flat prior on the model parameters, 9, so that the 
implied by the model, we have posterior is proportional to the likelihood in Equation (10): 


P(d\9) = n / dlij p(dij\Xij)p(\ij\9), (4) 

i>j J 

where 

(Vi jXi i) di -i 

P(dij\Xjj) = — — ^ exp (-VtjXij) (5) 

a i,j ■ 


is the Poisson probability of drawing dj j counts in a volume 
Vi j = i);j(100Mpc /i” 1 ) 3 when the underlying rate is Xjj per 
(lOOMpc/r 1 ) 3 . 

The distribution of the underlying rates implied by our model, 
p(Xj j \0). must be estimated from the observed ; . Applying 

Bayes’ rule again, we have 


p(hj\0) = P&ij\ n ij(0)) = 


P(njj\kij)p(hj) 

P( n i,j) 


(6) 


The counts observed in a model with underlying rate Xjj are 
Poisson distributed, so 


K! 


P(n uj \Xij) = — exp (-X U j). 


(7) 


We choose a Jeffreys prior 10 on Xjj, 


pO-i,j) = 


1 


whence 


P(hj\Q) = f- exp (-Xjj). 

r (riij + 5 ) 

Combining Equations (9) and (5), we find that 


p(d\9) = \ 


v tj T (s + di.j + >h,j) 

\,j (1 + Vij)i +di - i+ni ’ 1 dij !T ( j + «,j) 


( 8 ) 


(9) 


( 10 ) 


10 Note that the use of the Jeffreys prior implies that (A .,-j) = tiij + (1/2). A 
flat prior would have (A-jj) = nij + 1. Both of these priors produce 
well-defined likelihoods even when mj = 0 with dij 0. The 
maximum-likelihood estimator, p(A.,*j|0) = 8(kij — nij), while unbiased, 
produces likelihoods of zero if riij = 0, even if only a single count appears in 
that bin of the data (i.e., dij = 1). A prior that gives (A., j) = tiij (i.e., a prior 
that gives a distribution with unbiased mean) is = A./’j, which results in 

a non-normalizable likelihood when ntj = 0. 


p(9\d) cx p(d\9). (11) 

4. RESULTS AND DISCUSSION 

Figure 2 plots the likelihood values for each model used in 
this study and shows that likelihood values are very sensitive to 
model parameters and that only a few models are able to accu- 
rately reproduce the observed XLFs. Therefore, this comparison 
is very useful for eliminating regions of our model parameter 
space. Table 2 lists the 6 models with log(L(D|M)/L re f) > —3, 
where L ie f is the highest likelihood value among our 192 mod- 
els. These models are 205, 229, 277, 245, 253, and 273. These 
models all have low CE efficiencies (o'er = 0.1), a 50-50 mass 
ratio distribution, an IMF exponent of —2.7 (with the exception 
of model 273), and ?7 w i n d = 1. 0-2.0. Recall that our CE effi- 
ciency parameter really represents «(>; x /., so low values of 
o'er could alternatively be interpreted as these systems having 
a high envelope binding energy, which has been found to be 
true for massive stars (Dominik et al. 2012). The models that 
have 77 w i n d = 2.0, o'er. = 0.1, IMF exponent of —2.7, and a 
flat ^-distribution also have fairly high likelihood values (model 
numbers 13, 37, 61, and 85). 

It should be noted that our likelihood calculation takes 
into account the number of samples in each bin. The overall 
likelihood values are much more sensitive to bins with higher 
sample counts (i.e., those at lower luminosity). 

Figure 2 and Table 2 show that allowing/not allowing CE-HG 
phases has only a relatively small effect on the likelihoods of our 
best models. In addition, DC BH kicks only have an appreciable 
effect on likelihoods for models with lower wind mass-loss rates 
(»7wind = 0.25-1.0). For models with f; w j n d = 2.0, such as those 
that make up the majority of our top models, DC BH kicks 
have little effect. Thus, these two parameters are not very well 
constrained by our analysis. 

4.1. Comparison with F13 and T&G08 

F13 use the same PS models used in this work (with the 
inclusion of pure “twins” models) to study the evolution of 
the overall population of XRBs in the universe. They compare 
with X-ray observations of local galaxies that give estimates 
of the specific X-ray luminosity of XRBs in the local universe 
(Lehmer et al. 2010; Boroson et al. 2011; Mineo et al. 2012). 
They calculate the likelihood of each model based on these 
data and find the six highest likelihood models to be, in order 
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39 40 41 42 39 40 41 42 39 40 41 42 


log( L x (ergs s 1 ) [0.5-2. 0 keV] ) 



39 40 41 42 39 40 41 42 

log( L x (ergs s' 1 ) [0. 5-2.0 keV] ) 



39 40 41 42 39 40 41 42 

log( L x (ergs s' 1 ) [0. 5-2.0 keV] ) 

Figure 3. Cyan region shows the area bounded by the six highest likelihood models from this study, the gray checkered region shows that bounded by the six highest 
likelihood models from F13, the black points with error bars are data from T&G08, the solid black lines are from our highest likelihood model (205), and the dashed 
black lines are from the highest likelihood model from F13 (245). Top: total galaxy population. Bottom left: early-type galaxies. Bottom right: late-type galaxies. 
The XLFs from our highest likelihood models are very similar to those from the F13 models; however, they underproduce bright early-type galaxies and very bright 
(Lx > 10 41 erg s" ' 1 ) late-type galaxies compared with observations. Section 4.1 discusses the causes of these discrepancies in more detail. 

(A color version of this figure is available in the online journal.) 


of likelihood, 245, 229, 269, 205, 249, and 273 (see Table 2 
for model parameters and likelihood values from this study). 
Figure 3 compares the six highest likelihood models from 
this work with those of FI 3. Four out of our top six models 
are also among the six highest likelihood models from FI 3, 
so it is no surprise that the region bounded by the models 
in this work is very similar to that bounded by the models 
from F13. 

Figure 4 shows XLFs for our highest likelihood model 
(205), with and without hot gas emission, plotted against the 
data from T&G08. Figure 5 plots XLFs similar to Figure 4, 
but for the highest likelihood model from F13 (245). These 
plots show that our models are able to reproduce the redshift 
evolution of the observed XLFs. Consistent with the analysis of 
T&G08, the XLF evolution is driven almost entirely by late-type 
galaxies. Our models also reproduce the shape of the early-type 


XLF and the normalization of the late-type XLF, though they 
drastically underproduce bright early-type galaxies and they fail 
to reproduce the shape of the bright (L x > 10 41 erg s [ ) end of 
the late-type XLF. 

Figures 4 and 5 also show that hot gas can have a large effect 
on the shape and normalization of the XLF, showing that hot 
gas emission dominates the integrated X-ray luminosity of the 
brightest galaxies in our sample. Adding in the hot gas emission 
suppresses the redshift evolution for the early-type galaxy XLF. 
For galaxies with L x > 10 40 erg s 1 at low (z < 0.8) redshift, 
emission from XRBs accounts on average for only l%-5% and 
~15% of early- and late-type galaxy emission, respectively. 
However, as we will discuss in Section 4.3, the XLFs are still 
rather sensitive to changes in our model parameters, as seen in 
the varying likelihood values shown in Figure 2. The important 
role of hot gas emission on the XLF means that our simplistic 
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Figure 4. Luminosity function for the total galaxy population using our highest 
likelihood model, number 205. The solid lines show the XLFs from our models 
with both XRB and hot gas emission. The dot-dashed lines show XLFs for just 
XRB emission. The data points and associated error bars are taken from T&G08. 
Consistent with the analysis of T&G08, the overall XLF evolution is driven 
almost entirely by late-type galaxies. However, the model fails to reproduce the 
correct normalization of the early-type XLF and the shape of the observed XLF 
of very bright (Lx > 10 41 erg s _1 ) late-type galaxies. Section 4.1 discusses the 
possible causes for these discrepancies in more detail. Hot gas plays an important 
role in the normalization of the XLF, especially for early-type galaxies, where 
the hot gas contribution also affects the XLF evolution with redshift. 

(A color version of this figure is available in the online journal.) 


prescriptions for hot gas emission add a great deal of uncertainty 
to our models and could be a major source of the discrepancies 
between the models and observations, particularly for early- 
type galaxies. While our method is motivated by observations, 
it does not take into account the internal characteristics of the gas 
that contribute to its emission, such as density and metallicity. 
Further, the relations used for early-type galaxies were derived 
only from low-redshift sources, which may not be accurate for 
the high-redshift galaxies studied here. 

In addition, the Gil semi-analytic model underproduces 
massive galaxies at high redshift. This will lead to less large 
elliptical galaxies, which could explain part of our discrepancy 
at higher redshift. 

Another aspect of our models that can account for the 
underproduction of bright early-type galaxies is that our PS 



Figure 5. Same as Figure 4, but for model 245, the highest likelihood model 
from FI 3 and the fifth highest likelihood models from this work. These XLFs 
are similar to that of our highest likelihood model, 205. 

(A color version of this figure is available in the online journal.) 


models only take into account LMXBs formed in the field and 
not those formed dynamically in GCs. Dynamically formed 
LMXBs are believed to play a significant role in old, massive, 
GC-rich elliptical galaxies, (e.g., Humphrey & Buote 2008; 
Zhang et al. 2012). These LMXB populations can have a 
significant contribution to the integrated X-ray luminosity of 
bright early-type galaxies, as they can make up over half of the 
total number of LMXBs in a galaxy (Irwin 2005). So, including 
dynamically formed LMXBs in our models could increase 
the number of bright LMXBs, and therefore the total LMXB 
luminosity, in early-type galaxies by a factor of ~3. Changing 
the ^-distribution in model 245 from a 50-50 to a flat distribution 
is a good proxy for this effect because it increases the LMXB 
population without changing the distribution of their physical 
properties (see Section 4.3). F13 show that doing this increases 
the total luminosity from LMXBs at all redshifts by a factor 
of two. We find that changing to a flat ^-distribution increases 
the low-luminosity end of the early-type XLF by ~0.3 dex 
(see Figure 7). The effect of including dynamically formed 
LMXBs could have a similar but greater effect, bringing the low- 
luminosity end of the XLF closer to observations, but having 
little effect on higher luminosity galaxies. A more detailed 
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calculation will require information on the GC population of 
each galaxy, which is not included in the G1 1 catalog. 

For younger, star-forming galaxies, our models also have only 
a very basic formula to simulate starburst activity, which can 
occur, e.g., due to galaxy mergers. This would have a significant 
effect on the HMXB populations present in late-type galaxies, 
and the effect would not necessarily be constant with redshift, 
as merger rates may evolve with time (e.g., Conselice 2006). 
Thus, a more detailed SFH is needed to more accurately model 
HMXB populations of late-type galaxies. 

In addition, the higher end of the observed late-type galaxy 
XLF is more at risk from AGN contamination, even with the 
efforts of T&G08 to minimize this effect. Since the observations 
of T&G08, the depth of the X-ray data, combined with better 
multiwavelength data, has allowed for more accurate classifi- 
cations of the X-ray sources. Of the 56 1 Ms CDF-S sources 
used in T&G08, we find 53 counterparts with 4 Ms exposure 
using a matching radius of 2.5 arcsec. The missing three sources 
may have been false-positive sources in the 1 Ms data. Of these 
53 sources, we find that 25 of them are classified as normal 
galaxies and 28 of them as AGN according to the 6 criteria 
highlighted in Section 3.1 of Lehmer et al. (2012). Therefore, 
it is possible that the T&G08 data points will be lowered by 
~0.3 dex. However, it is difficult to know in detail how this 
affects the TG08 luminosity functions and recomputing the lu- 
minosity functions is beyond the scope of this work. 

4.2. High-redshift Predictions 

Figure 6 plots the X-ray luminosity density from normal 
galaxies as derived from our highest likelihood model, 205. 
The overall evolution (black line in Figure 6) is very similar to 
that of the observed SFH of the universe. It is also similar to the 
evolution of the specific XRB X-ray luminosity of the universe 
predicted in F13, despite the inclusion here of hot gas emission. 
This is evidence that XRBs drive the overall evolution of the 
normal galaxy X-ray luminosities out to at least z = 4. However, 
our predicted X-ray luminosity density reaches a maximum at 
z ~ 2.5, which is lower compared with the XRB models in 
F13 that reach a maximum at z ~ 3. This can be attributed 
to the inclusion of hot gas emission in our models, which has 
already been shown to have a noticeable effect on the shape and 
evolution of our XLFs, particularly for early-type galaxies. 

Splitting the galaxies into three luminosity bins, we find that 
the evolution of low (10 39 < Lx < 10 40 ergs -1 ) luminos- 
ity galaxy emission is small compared with the evolution for 
higher luminosity galaxies, which varies by an order of mag- 
nitude on the range of z = 0-4. The most luminous galax- 
ies (Lx > 10 41 ergs -1 ) reach a maximum around z = 3, 
which also approximately corresponds to the time of maxi- 
mum SFR density in the universe. Galaxies in the range of 
10 4 ° < Lx < 10 41 erg s -1 reach a maximum around z = 2. 
In the local universe, the low-luminosity galaxies dominate the 
normal galaxy X-ray emission. We do not go to higher red- 
shift here because our hot gas emission prescription relies on 
galaxy morphology, which becomes harder to classify at higher 
redshifts (van den Bergh 2002). 

4.3. Effects of Parameters on XLFs 

Figure 7 shows XLFs for different models compared to model 
245. Each model is chosen to encapsulate the effect that each 
parameter has on the shape of the XLF. Several parameters have 
significant effects on the shape of the XLFs, while others have 



Figure 6. X-ray luminosity density of normal galaxies vs. redshift for model 
205. The back line shows the contribution from all galaxies and the black 
squares show the total X-ray luminosity density from observations presented 
in T&G08. The X-ray luminosity density derived from our models is a factor 
of 2-3 lower than the observations. This is due to the fact that our models, 
relative to observations, underproduce bright early-type galaxies and very bright 
(Lx > 10 41 ) late-type galaxies. Our model follows a similar evolution as the 
SFR density and the X-ray luminosity density from XRBs predicted in FI 3, 
indicating that XRBs continue to drive the evolution of the normal galaxy X-ray 
emission at higher redshifts. However, our models reach a maximum luminosity 
density at z ~ 2.5, compared with z ~ 3 for the XRB models in FI 3. This can 
be attributed to the inclusion of hot gas emission, which has been shown to have 
an appreciable effect the evolution of our XLFs when compared with XRB only 
models. The red, green, and blue lines plot the specific X-ray luminosity for 
different luminosity bins. The amount X-ray emission from lower luminosity 
galaxies evolves much less with redshift than that from brighter galaxies. In the 
local universe, most of the normal galaxy X-ray emission comes from lower 
luminosity galaxies. 

(A color version of this figure is available in the online journal.) 

only minimal effects. Model 245 was chosen because it has 
both a high likelihood and is more sensitive to certain changes 
in parameters, thus better illustrating the different effects on 
our XLFs. 

The CE efficiency parameter (q!ce) dictates how efficiently 
orbital energy is converted to thermal energy that will expel the 
envelope. A lower efficiency means that it will take more orbital 
energy to expel the envelope. This parameter mainly affects 
LMXBs, as most LMXBs formed in the field must go through 
a CE phase. The CE phase plays an important role in making 
the orbit close enough to allow for Roche-lobe overflow (RLO), 
but a lower CE efficiency leads to even more orbital decay and 
a higher rate of mergers, overall decreasing the rate of LMXB 
formation. This effect can be seen by comparing models 245 
and 248. HMXBs are not as strongly affected by changes in 
o'er, as they have other formation channels available that do not 
involve CE phases (Linden et al. 2010; Valsecchi et al. 2010). 
This parameter mainly affects the lower luminosity end of the 
XLF, where a higher CE efficiency increases the number of 
bright galaxies due to an increased LMXB population. 

Wind mass-loss rates affect the evolution of high-mass stars 
in two major competing ways. Higher wind mass-loss rates will 
increase the accretion rates of wind-fed HMXBs, increasing 
their luminosity. On the contrary, lower winds will result 
in a lower overall mass loss of the primary star and hence 
increase the formation rate of massive BHs. BH-XRBs tend 
to be more luminous than XRBs with NS accretors. This 
is because, on one hand, they can form stable RLO XRB 
systems with massive companions and, on the other hand, 
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Figure 7. Highest likelihood model from FI 3 (245, black line) compared with 
other models to illustrate the effects of different parameters on the shape of the 
XLF. Different parameters have varying effects on the shape of the XLFs. 
Shown here are models with higher acE (248), CE-HG allowed (269), no 
BH kicks (197), steeper IMF (241), lower wind mass loss (261), and a flat 
distribution (53). 

(A color version of this figure is available in the online journal.) 


BHs drive higher accretion rates due to their higher mass and, 
therefore, higher Eddington limits. In this way, weaker stellar 
winds can increase the luminosities of both LMXB and HMXB 
populations. Comparing models 245 and 26 1 , we see that weaker 
stellar winds increase the number of bright early- and late-type 
galaxies, so the latter effect is dominant. While lower stellar 
winds will help our highest likelihood model match observations 
of early-type galaxies, it would overproduce bright late-type 
galaxies. 

Changing the initial binary mass ratio between a flat distri- 
bution and a 50% twins and 50% flat distribution affects both 
the HMXB and LMXB populations. As stated earlier, the binary 
mass ratio affects the secondary star in the system, which will 
eventually become the donor star in most cases. All XRBs are 
accreting mass onto a compact object, which can only come 
from a high-mass progenitor. Most LMXBs require a high ini- 
tial mass ratio in order to ensure a high-mass primary star that 
will evolve into a BH or NS with a lower mass companion. The 
“50-50” distribution adopted here forces mass ratios close to 
1 for half the binary population. This will decrease the LMXB 


population while increasing the HMXB population, as HMXBs 
require ratios closer to 1. Changing from “50-50” to a flat dis- 
tribution (comparing models 245 and 53) increases the lower 
luminosity end of the early- and late-type galaxy XLLs and 
slightly decreases the number of very bright late-type galaxies 
galaxies. 

Allowing small natal kicks for DC BHs affects the HMXB 
and LMXB populations in two competing ways. On one hand, 
natal kicks can enhance the formation of RLO-HMXBs (Linden 
et al. 2010). On the other hand, natal kicks inject energy 
into the binary system and could result in the widening or 
the complete disruption of the system, thereby decreasing the 
formation of HMXB and LMXB with BH accretors. Linden 
et al. (2009) find that imparting small natal kicks to DC BHs is 
necessary in order to reproduce the lack of observed wide orbit 
BH-XRBs. Comparing models 245 and 197, not allowing natal 
kicks increases the high-luminosity end of the late-type galaxy 
XLL and has little effect on the early-type galaxy population. 

Within our grid of PS models, we also have the IMP power- 
law exponent as a free parameter, allowing it to be either —2.35 
or —2.7. It is instructive to note that the IMP referred to in this 
work represents the integrated galaxy IMP (Weidner & Kroupa 
2005; Kroupa & Weidner 2003) and that this work only probes 
the high-mass region of the IMP, since the primary stars that are 
created via sampling the IMP must be massive enough to form a 
BH or NS. The slope of the power law at the high-mass end of the 
IMP affects the population of XRBs in a way similar to stellar 
winds, in the sense that a flatter IMP will produce relatively more 
massive BHs compared to a steeper one. A flatter IMF will result 
in more bright LMXBs and HMXBs. This effect can be seen 
by comparing models 245 and 241. As expected, a flatter IMF 
results in a higher number of bright galaxies, though the very 
high luminosity end of the early-type galaxy XLF is not strongly 
affected by this parameter. This is not surprising, as Figure 4 
shows that hot gas emission dominates the high-luminosity end 
of the XLF for early-type galaxies. 

Finally, we found that allowing or not for all possible 
outcomes in CE phases with donor stars in the HG has only a 
small effect on the galaxy XLFs, slightly increasing the number 
of bright late-type galaxies (compare models 245 and 269). 
Thus, although this parameter affects the shape of the XLF of 
individual XRBs in a galaxy (Luo et al. 2012), it has a negligible 
effect on the integrated X-ray luminosity of a galaxy. 

The models that we find to agree best with observations have 
DC BH natal kicks, low CE efficiency, steep IMFs, ?? w i n d = 
1 .0-2.0, and a 50-50 mass ratio distribution. As outlined above, 
these parameters make for a limited population of both LMXBs 
and BH-XRBs. The HMXB population is limited by the steeper 
IMF, BH natal kicks, and higher winds, but also strengthened 
by the 50-50 mass ratio distribution. However, the latter effect 
on the XLF is much weaker, as shown in Figure 7. 

5. CONCLUSIONS 

Using data from the Millennium Cosmological Simulation 
and the semi-analytical analysis conducted by Gil in tandem 
with the binary PS code, StarTrack, we simulated the population 
of XRB s within normal galaxies in a large volume of the universe 
from z — 0 to ~20. Assuming that galaxy X-ray emission is 
solely due to XRBs and hot gas, we calculated the integrated 
X-ray luminosity of each galaxy in this cosmic volume and 
compared the resulting galaxy XLFs to the observational XLFs 
of T&G08. 
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In this paper, we presented data from 192 binary PS models, 
varying parameters that have the largest effect on binary star 
evolution (see Table 1 for a list of the parameters). We use 
a likelihood calculation method to compare each model with 
the results from T&G08. From this analysis we find that our 
theoretical XLFs are sensitive to many of our model parameters 
and that only a few of our models are able to reproduce the 
most recent observations of X-ray bright normal galaxies. Our 
highest likelihood models are also among the highest likelihood 
models from a separate analysis presented in FI 3. This confirms 
that our results are consistent with their separate analysis, which 
compares these same models with the observed overall emission 
from LMXB and HMXB populations in the local universe. To 
have only ~10 models from our 192 model grid best match 
observations in 2 separate analysis shows that we are able to 
provide self-consistent constraints on the XRB parameter space. 

We find that our highest likelihood models are those with 
a lower LMXB population due to a low CE efficiency and a 
50-50 mass ratio distribution, and a lower BH-XRB and HMXB 
population due to higher winds, and a steeper IMF. Our models 
do well in reproducing the normalization and evolution of the 
total and late-type galaxy XLFs, as well as the evolution and 
shape of the early-type XLF. 

Our models show that hot gas emission has a large effect on 
the shape of the XLFs, and it significantly affects the redshift 
evolution of the early-type galaxy XLF, causing it to remain 
nearly constant out to z = 1.4. 

We show that the observed redshift evolution of the normal 
galaxy XLF continues out to higher redshift, with the specific 
normal galaxy X-ray luminosity evolving in a way similar to the 
SFH of the universe and consistent with the evolution of XRB 
emission found in F13. This is evidence that the XLF evolution 
is driven by XRB evolution even out to higher redshifts. Our 
models also show that hot gas emission causes the point of 
maximum normal galaxy X-ray luminosity density to shift to 
lower redshift compared with the XRB models in FI 3. 

However, despite these many successes, our models do not 
perfectly reproduce the observed XLFs. In particular, they fail to 
reproduce the observed normalization of the early-type galaxy 
XLF, greatly underestimating the number of bright early-type 
galaxies. Our highest likelihood models also fail to reproduce 
the shape of the high (L x > 10 41 erg s') luminosity end of the 
late-type galaxy XLF, particularly for higher redshifts. 

Our models have limitations that may have caused these dis- 
crepancies. For one, we do not take into account dynamically 
formed LMXBs, which could significantly increase the nor- 
malization of the model early-type galaxy XLF. For late-type 
galaxies, the XRB luminosities have a higher contribution from 
HMXB populations, which are very sensitive to evolving SFHs. 
However, the SFHs used from the G1 1 catalog are limited in 
their detail and our method for simulating the effect of star- 
bursts is very rudimentary. In addition, our prescription for hot 
gas, though based on observations, is very basic and could add 
inaccuracy to our X-ray luminosities as well as the selection 
of our best-fitting models. A more detailed model is needed 
to more accurately model the hot gas emission, particularly in 
early-type galaxies. 

In addition to limitations in our models, the observations 
of very bright galaxies are subject to the possibility of AGN 
contamination, which could artificially increase the observed 
high-luminosity data points from T&G08. 

Despite these shortcomings, this work represents a first 
careful attempt to study how XRBs control the L x distributions 


of different types of galaxies. As such, it provides an important 
theoretical base for future X-ray observations of normal galaxies 
at high redshift. It also shows that XRB populations are closely 
linked with the growth of galaxies. This work lays the ground 
for future work using X-ray observations and cosmological 
simulations of galaxies to provide a new way to constrain our 
models of binary evolution, as well as study the role played 
by XRBs in galaxy formation and evolution through feedback 
processes. 
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